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The recently completed calculation of the full electroweak 0(ce) corrections to the charged-current four-fermion pro- 
duction processes e+e — — > v t t+ fi~ 9^, ud/x — P M , and udsc is briefly reviewed. The calculation is performed using 
complex gauge-boson masses, supplemented by complex couplings to restore gauge invariance. The evaluation of the 
occurring one-loop tensor integrals, which include 5- and 6-point functions, requires new techniques. The effects of the 
complete 0(a) corrections to the total cross section and to some differential cross sections of physical interest are dis- 
cussed and compared to predictions based on the double-pole approximation, revealing that the latter approximation 
is not sufficient to fully exploit the potential of a future linear collider in an analysis of W-boson pairs at high energies. 



1. INTRODUCTION 



At LEP2, W-pair-mediated four-fermion (4/) production was experimentally explored with quite high precision 
(see Ref. and references therein). The LEP2 measurements had set the scale in accuracy in the theoretical 
predictions for W-pair-mediated 4/ production, as reviewed in Refs. 0,0]. For LEP2 accuracy, it was sufficient to 
include corrections in the so-called double-pole approximation (DPA), where only the leading term in an expansion 
about the poles in the two W-boson propagators is taken into account. Different versions of such a DPA have been 
used in the literature 0, IE 0, 0, 0] • Although several Monte Carlo programs exist that include universal corrections, 
only two event generators, YFSWW 0,0] and RacoonWW 0,0 Q, include non-universal corrections. 

In the DPA approach, the W-pair cross section can be predicted within ~ 0.5% (0.7%) in the energy range 
between 180 GeV (170 GeV) and ~ 500 GeV, which was sufficient for the LEP2 accuracy of ~ 1%. In the threshold 
region (y/s < 170 GeV), the DPA is not reliable, and the best available prediction results from an improved Born 
approximation (IBA) based on leading universal corrections only, and thus possesses an intrinsic uncertainty of ~ 2%. 

At a future International e + e _ Linear Collider (ILC) EH HQ, 

the accuracy of the cross-section measurement 
will be at the per-mille level, and the precision of the W-mass determination is expected to be ~ 10 Me V I l4l by 
direct reconstruction and ~ 7 MeV from a threshold scan of the total W-pair-production cross section [n[ Il2j. The 
theoretical uncertainty (TU) for the direct mass reconstruction at LEP2 is estimated to be of the order of ~ 5 MeV 
to < 10 MeV 0, based on results of YFSWW and RacoonWW; theoretical improvements are, thus, desirable 
for an ILC. For the cross-section prediction at threshold the TU is ~ 2%, because it is based on an IBA, and thus is 
definitely insufficient for the planned precision measurement of Mw in a threshold scan. The TU in constraining the 
anomalous triple-gauge-boson coupling A 7 was estimated to be ~ 0.005 0] for the LEP2 analysis. Since a future 
ILC is more sensitive to anomalous gauge-boson couplings than LEP2 by more than an order of magnitude, a further 
reduction of the uncertainties resulting from missing radiative corrections is necessary. 

Recently we have completed the first O(a) calculation (improved by higher-order ISR) for the 4/ final states 
v t t + fi~ v^, vAu,~u,i,, and udsc, which are relevant for W-pair production. We have presented results on total cross 
sections in Ref. 01 an d on various differential distributions in Ref. 2Q] . The latter publication also contains technical 



PSN 0609 



details of the actual calculation, which is rather complicated. 1 In the following we briefly describe the salient features 
of the calculation and present a selection of numerical results that are relevant for a future ILC, comprising total 
cross sections and some phenomenologically interesting distributions. 

2. METHOD OF CALCULATION 

The actual calculation builds upon the RacoonWW approach [7j, where real-photonic corrections are based on 
full matrix elements and virtual corrections are treated in DPA. Real and virtual corrections are combined either 
using two-cutoff phase-space slicing or employing the dipole subtraction method |2l| for photon radiation. We also 
include leading- logarithmic ISR beyond O{o) in the structure-function approach (Ref. Q and references therein). 

2.1. Technical issues 

In contrast to the DPA approach, the one-loop calculation of an e + e~ — > 4/ process requires the evaluation of 5- 
and 6-point one-loop tensor integrals. For the generic /1/2/3/4 final state, where /1 and are different fermions 
excluding electrons and electron neutrinos and ]i and f± their isospin partners, there are 40 hexagon diagrams, 
112 pentagon diagrams, and 227 (220) box diagrams in the conventional 't Hooft-Feynman gauge (background-field 
gauge). A survey of Feynman diagrams can be found in Ref. J20j. We calculate the 6-point integrals by directly 
reducing them to six 5-point functions, as described in Refs. [i2,|2g. The 5-point integrals are reduced to five 4-point 
functions following the method of Ref. [3 . Note that this reduction of 5- and 6-point integrals to 4-point integrals 
does not involve inverse Gram determinants composed of external momenta, which naturally occur in the Passarino- 
Vcltman reduction 25] of tensor to scalar integrals. The latter procedure leads to serious numerical problems when 
the Gram determinants become small, which happens usually near the boundary of phase space but can also occur 
within phase space because of the indefinite Minkowski metric. 

We use, however, Passarino-Veltman reduction to calculate tensor integrals up to 4-point functions, which involves 
inverse Gram determinants built from two or three momenta. This, in fact, leads to numerical instabilities in phase- 
space regions where these Gram determinants become small. For these regions we have worked out two "rescue 
systems": one makes use of expansions of the tensor coefficients about the limit of vanishing Gram determinants; 
in the other, alternative method we numerically evaluate a specific tensor coefficient, the integrand of which is 
logarithmic in Feynman parametrization, and derive the remaining coefficients as well as the scalar integral from it 
algebraically. This reduction does not involve inverse Gram determinants. 

In addition to the evaluation of the one-loop integrals, also the evaluation of the three spinor chains corresponding 
to the three external fermion-antifermion pairs is non-trivial, because the chains are contracted with each other 
and/or with four-momenta in many different ways. There are O(10 3 ) different chains to calculate, so that an 
algebraic reduction to a standard form which involves only very few standard chains is desirable. We have worked 
out algorithms that reduce all occurring spinor chains to O(10) standard structures without introducing coefficients 
that lead to numerical problems. These algorithms are described in Ref. in detail. 

2.2. Conceptual issues 

The description of resonances in (standard) perturbation theory requires a Dyson summation of self-energy inser- 
tions in the resonant propagator in order to introduce the imaginary part provided by the finite decay width into the 
propagator denominator. This procedure in general violates gauge invariance, i.e. destroys Slavnov-Taylor or Ward 



1 Recently the authors of the GRACE/1-LOOP system reported on progress towards a full one-loop calculation for e+c — » fi ;/ M ud 
in Ref. Ilgl so that one can expect that the system will be able to deal with e + e _ — > 4/ processes at one loop in the near future. 



PSN 0609 



identities and disturbs the cancellation of gauge-parameter dependences, because different perturbative orders are 
mixed (see, for instance, Ref. Q and references therein). The DPA provides a gauge-invariant answer in terms of an 
expansion about the resonance, but in the full calculation we are after a unified description that is valid both for reso- 
nant and non-resonant regions in phase space, without any matching between different treatments for different regions. 

For our calculation we have generalized |2Q| the so-called "complex-mass scheme" , which was introduced in Ref. |9J 
for lowest-order calculations, to the one-loop level. In this approach the W- and Z-boson masses are consistently 
considered as complex quantities, defined as the locations of the poles in the complex p 2 plane of the corresponding 
propagators with momentum p. Gauge invariance is preserved if the complex masses are introduced everywhere in 
the Feynman rules, in particular, in the definition of the weak mixing angle, which is now derived from the ratio of 
the complex masses. The (algebraic) relations, such as Ward identities, that follow from gauge invariance remain 
valid, because the gauge-boson masses are modified only by an analytic continuation. Since this continuation already 
modifies the lowest-order predictions by changing the gauge-boson masses, double-counting of higher-order effects 
has to be carefully avoided by an appropriate renormalization procedure. 

The use of complex gauge-boson masses necessitates the consistent use of these complex masses also in loop 
integrals. To this end, we have derived all relevant one-loop integrals with complex internal masses. The IR-singular 
integrals were taken from Ref. j^. Concerning the non-IR singular cases, we have analytically continued the results 
of Ref. for the 2-point and 3-point functions, and the relevant results of Ref. [3] for the 4-point functions. We 
have checked all these results by independent direct calculation of the Feynman-parameter integrals. 



2.3. Checks on the calculation 

In order to prove the reliability of our results, we have carried out a number of checks, as described in more detail 
in Ref. ^^|. We have checked the structure of the (UV, soft, and collinear) singularities, the matching between 
virtual and real corrections, and the gauge independence (by performing the calculation in 't Hooft-Feynman gauge 
and in the background-field gauge |29|). The most convincing check for ourselves is the fact that we worked out the 
whole calculation in two independent ways, resulting in two independent computer codes the results of which are in 
good agreement. All algebraic manipulations, including the generation of Feynman diagrams, have been done using 
independent programs. The amplitudes are generated with FeynArts, using the two independent versions 1 and 3, 
as described in Refs. and [il]], respectively. The algebraic manipulations are performed using two independent 
in-house programs implemented in Mathematica, one of which builds upon FormCalc 0] ■ For the calculation 
of the loop integrals we use the two independent in-house libraries which employ the different calculational methods 
sketched above for the numerical stabilization. 



3. NUMERICAL RESULTS 

The precisely defined input for the numerical results presented in the following can be found in Refs. jlflL bflj]. 

Figure^Jdepicts the total cross section for the energy ranges of LEP2 and of the high-energy phase of a future ILC, 
focusing on the leptonic final state v t t + D^. The respective figures for the relative corrections 5 to the semileptonic 
and hadronic final states look almost identical, up to an offset resulting from the QCD corrections. Specifically, the 
upper plots show the absolute prediction for the cross section including the full 0(a) corrections and improvements 
from higher-order ISR. The lower plots compare the relative corrections as obtained from the full 0(a) calculation, 
from an IBA, and from the DPA. The IBA 10] implemented in RacoonWW is based on universal corrections only 
and includes solely the contributions of the CC03 diagrams. The DPA of RacoonWW comprises also non-universal 
corrections |j| and goes beyond a pure pole approximation in two respects. The real-photonic corrections are based 
on the full e + e~ — > 4/ + 7 matrix elements, and the Coulomb singularity is included for off-shell W bosons. Further 
details can be found in Ref. Q- 

A comparison between the DPA and the predictions based on the full 0(a) corrections reveals differences in the 
relative corrections S of < 0.5% (0.7%) for CM energies ranging from y/a ~ 170 GeV to 300 GeV (500 GcV). This 
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Figure 1: Absolute cross section a (upper plots) and relative corrections S (lower plots) to the total cross section without cuts 
for e + e~ — ► ^ r r + /i _ i/ M obtained from the IBA, DPA, and the full O(a) calculation (ee4f). All predictions are improved by 
higher-order ISR. (Taken from Ref. Q.) 



is in agreement with the expected reliability of the DPA, as discussed in Refs. |ij |6j, |7j . At higher energies, the 
deviations increase and reach 1—2% at \/s = 1— 2TcV. In the threshold region (y/s < 170 GeV), as expected, the 
DPA also becomes worse w.r.t. the full one-loop calculation, because the naive error estimate of (a/ir) x (T\y/M\y) 
times some numerical safety factor of 0(1—10) for the corrections missing in the DPA has to be replaced by 
(a/ir) x Tyj/(y/s — 2Mw) in the threshold region and thus becomes large. In view of that, the DPA is even 
surprisingly good near threshold. For CM energies below 170 GeV the LEP2 cross section analysis was based 
on approximations like the shown IBA, which follows the full one-loop corrections even below the threshold at 
s/s = 2Mw within an accuracy of about 2%, as expected in Ref. ^jj]. More results on total cross sections, including 
numbers on leptonic, semileptonic, and hadronic final states, can be found in Ref. 0]. 

The distributions in the invariant mass of the W + boson and in the cosine of the W + production angle 6>w+ are 
shown in Figure [21 for the process e + c~ — * ud/i~P M at s/s = 500 GeV. Further distributions, also for y/s — 200 GeV, 
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Figure 2: Distributions in the invariant mass of the W + boson (l.h.s.) and in the cosine of the W + production angle with 
respect to the e + beam (r.h.s.) and the corresponding corrections (lower row) at s/s = 500 GeV for e + e" -> uduji^. The 
inset plot shows the relative difference A between the full O(a) corrections and those in DPA. (Taken from Ref. (2£j] .) 



are presented in Ref. For the invariant-mass distribution the full O(o) calculation and the DPA agree within 
~ 1% both for LEP2 and ILC energies. For the W-production-angle distribution this is also the case in the LEP2 
range (see Fig. 12 of Ref. but at 500 GeV the difference of the corrections in DPA and the complete 0(a) 

corrections rises from —1% to about —2.5% with increasing scattering angle. Note that such a distortion of the shape 
of the angular distribution can be a signal for anomalous triple gauge-boson couplings. 
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